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TIRE PRESSURE LOSS DETECTION 
BACKGROUND 

[0001] The present invention relates to air pressure loss detection in one or 

more tires of a vehicle. 

[0002] During the operation of a vehicle, significant loss of pressure in one or 

more tires may cause the driver to lose control of the vehicle or a tire failure. 
Moreover, event slight pressure losses in a tire can result in reduced service life for 
that tire. Therefore, it is desirable to be able to inform the driver that one or more 
tires are losing pressure, especially before hazardous conditions are reached. 
[0003] Certain systems measure the pressure of each tire directly and relay 

this information to the driver. Other systems measure changes in the effective 
rolling diameter of the tires caused by losses in pressure or its response to road 
induced vibration. Yet in other systems the rotational velocities of the tires are 
measured under certain operating situations to identify conditions of the tires. For 
instance, loss of air in a tire may cause an increase in rotational resistance of the 
tire, and hence a perceptible change in the tire rotational velocity. 
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BRIEF SUMMARY OF THE INVENTION 

[0004] As a vehicle maneuvers over a road surface, ground excitation 
generates angular speed variations or torsional vibration of one or more wheels of 
the vehicle caused by for-aft displacements of the wheels. The frequency of the 
speed variations depends on the tire sidewall stiffness, which is mainly influenced by 
the tire pressure. The angular speed variation of a wheel shifts from a higher 
frequency to a lower frequency at a given speed when there is tire pressure loss. 
[0005] In accordance with the invention, the loss of pressure in one or more 

tires is monitored by detecting the angular speed variations of one or more wheels of 
the vehicle and analyzing the frequency of the speed variations over a specified 
wheel revolution period. Changes in the frequency are related to pressure loss in 
one or more tires, which is indicated to the driver, for example, by displaying the tire 
pressure loss information on a display. 

[0006] Other features and advantages will be apparent from the following 

drawings, detailed description, and claims. 

BRIEF DESCRIPTION OF THE DRAWINGS 

[0007] FIG. 1 illustrates an ABS encoder employed to monitor shifts in 

torsional vibration of an axle in accordance with an embodiment of the invention. 
[0008] FIG. 2 depicts the frequency spectrum for for-aft torsional vibrations of 

an axle. 

[0009] FIG. 3 depicts the distortion of a frequency spectrum caused by pole 

pitch errors. 

2 



Attorney Docket No. 10543-074 

[0010] FIG. 4 depicts a frequency spectrum for a tire pressure loss of about 

1.2 bars after eliminating the pole pitch error for a vehicle driven on a regular road. 
[0011] FIG. 5 depicts a frequency spectrum for a tire pressure loss of about 

1.2 bars after eliminating the pole pitch error for a vehicle driven on a rough road. 
[0012] FIG. 6 depicts frequency spectrum of data for 1.5 bars and 2.7 bars 

and respective curve fitted results. 

[0013] FIG. 7 depicts the frequency spectra of a front left tire for various 

pressure losses. 

[0014] FIG. 8 depicts the frequency spectra of a right left tire for various 

pressure losses. 

[0015] FIG. 9 depicts the frequency spectra of a front right tire for various 

pressure losses. 

[0016] FIG. 10 depicts the frequency spectra of a rear right tire for various 

pressure losses. 

[0017] FIG. 11 depicts a state diagram for torsional vibration detection in 

accordance with an embodiment of the invention. 

[0018] FIG. 12 depicts a flow diagram of a sequence of steps for the torsional 

vibration detection in accordance with an embodiment of the invention. 
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DETAILED DESCRIPTION 

[0019] The following terminology is helpful for understanding various aspects 
of the present invention: 
FFT - Fast Fourier Transform 

DFFT - Discrete Finite Fourier Transform 

© - A finite period of angle (as a multiple of 2n radian ) for estimating the Fourier 
transform 

TV - A finite number of sampled points for estimating the DFFT 
0 n - Angle of the n th signal sampled 

AO - Angle interval between two consecutive samples, 

(i.e. A^=— {In radian ) ) 
TV 

At n - A period of time during the n ih rotation of AO 

r{6) - A signal sampled in the angle 0 -nAO (i.e. the time period between two 
encoder flanks of rising edges or down edges; r{0)=r{nA0)=At n ) 

M - Number of frequencies detected in a frequency spectrum 

AQ - Angular frequency resolution in the angular domain that is equal to the 
difference between consecutive frequencies {cycles I 2k radian) 

Q - Angular frequency {cycles tin radian) defined in the angular spectrum 

Q = mAQ w = 0,l ? 2 ? 3,... 5 M-l 
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co - Angular velocity of an axle ( 2 n radian I sec ) 

a / - Frequency resolution in the time domain or difference between consecutive 
frequencies ( Hz ) 

[0020] During the operation of a vehicle, for example, ground excitation 

generates torsional vibration or angular speed variations of one or more wheels or 
other powertrain components of the vehicle, such as the axles of the driven wheels. 
For example, in FIG. 2 illustrates the frequency spectrum of the vibration caused by 
such excitations as a vehicle moves along a straight line at speeds at or above 40 
kph. As shown, the torsional vibration caused by for-aft excitation in the direction of 
travel of the vehicle is about 40 Hz. The frequency of the torsional vibration 
depends on the tire sidewall stiffness, which is mainly influenced by the tire 
pressure. The torsional vibration of a wheel shifts from a higher frequency to a lower 
frequency at a given speed when there is a pressure loss in a tire. 
[0021] In accordance with various embodiments of the invention, a process or 

algorithm (see, e.g., FIGs. 11 or 12) detects the tire air state based on the extracted 
natural frequency of the torsional vibration or angular speed variations of the wheel 
on which the tire is mounted. The process employs DFFT analysis for analyzing the 
torsional vibration. A DFFT analysis is disclosed in U.S. Patent No. 6,374,163, the 
contents of which is incorporated herein in its entirety. Accordingly, in the present 
invention, vibration analysis through ABS wheel speed sensors enables detecting 
the torsional frequency shifting and hence tire pressure loss. 

[0022] Referring now to FIG. 1, there is shown an ABS encoder 10 with N 

teeth 12, where during the rotation of the ABS encoder A/ f and A t r represent 



5 



Attorney Docket No. 10543-074 

the time period between two encoder flanks of rising edges 14 and down edges 16, 
respectfully. As the wheel rotates, appropriate software, implemented, for example, 
in an electronic control unit (ECU) of the vehicle, collects the number of tics passing 

one encoder flank (one tooth). Discontinuities in A / f or At r signals may be 

caused by speed changes, production/mounting variations of encoder, and/or rough 
roads. The ABS encoder (or tone wheel) produces multiple peak outputs per wheel 
revolution and is typically based on variable reluctance or Hall effect devices. The 
ABS encoder is provided to accurately measure wheel speeds for use as inputs for 
variable anti-lock brakes (ABS), dynamic stability program (DSP) control (and other 
yaw stability and roll-over stability systems), and traction control systems (TCS). 
[0023] A function is sampled at every spaced interval angle and the Fourier 
Transform is applied in the angular domain s(Q) , which is a function of the inverse 
angle e (i.e., q = l / e ) expressed in cycles per 2;rradian (i. e. cycles! '2 n), as 
given by the expression 

s(Q) = I" r{6) e- J27rno d0 (1) 

[0024] The Fourier integral transforms can be rewritten for discrete sample 
data systems. In the angular domain, 0 is the sampling angle and N is the number 
of samples taken during a period of 0 . In general, 0 is chosen as a multiple of in 
radian. Accordingly, the spectrum is given by the expression: 

s(n?AQ ) = — V r(nA 6 ) e - J2 * mn/N m = 0,1,2,..., — (2) 
N f? 0 2 

where the angle interval AO is considered a constant value, namely, 
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In 

A9 = — (3) 

Number of Teeth ABS encoder 



and the rotation angle is 

6 = n A 0 (4) 
The angular frequency resolution is provided by the expression 

AQ = 2 - cycles! 2 n (5) 
0 

where the angular frequency q is represented as w ao . 

Accordingly, the Fourier integral transforms a discrete sampled angular domain to 
the angular frequency domain in units of cycles! in. 

[0025] The spaced interval in time At n represents a period of time during the 

n th rotation of Ad , and At n in general is not a constant because of the variations of 

the angular velocity of a given wheel. 

[0026] If the sequence of sampled values is 

r(nA0) = At n (6) 

then from equation (2), the DFFT algorithm uses the following expression which 
applies equally to series of angles in real space: 

s(mACl ) = X At n e -i 2 * mnlN m = 0,1,2,..., ^ (7) 
Here, a contribution of spaced interval in time At n is At n e -J 2 * m » tN . 
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[0027] For the analysis of the torsional vibration, the sampling length N is 

calculated by the equation: 

[0028] A continuous contribution of each encoder flank 

At n (n = 0,1,2,3, ...,7V-1) to all the predefined frequencies is 
At n x e -j**»»'" m = 0,1,2,..., m - i. Note that, in general, more than one 
encoder flank is rotated during a software loop, and all the contributions of these 
encoder flanks to the spectrum are completed within the same software loop. 
[0029] In equation (2) the task of calculating the frequency spectrum is 

divided into N sub-tasks a/„ xe - J27Tmn/N n = 0,1,2,..., n - 1 . Hence, the total 

contribution for At n (n = 0,1,2,3,.. jV-1) creates a single angular frequency spectrum: 

0 N --> 



s 



( m AQ ) = — £ A t n x m = lj2) 3, M - 1 



The total tics ^ At n are accumulated to calculate an angular velocity co . 



n=0 



[0030] Rewriting equation (7) in complex space, the frequency spectrum is 

given by: 



s (m AQ) = — 

N 



f N-\ 



2n mn\ 



n C0S — N " J ^ " 



sin 



J n = 0 



2 n m n 
~N 



(9) 



n i n N 
m = 0,1,2, ... ,— 



and the amplitude of the angular frequency is: 
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■ , ™. 0 a (2xmn\) 2 ffc' . . ( 2 * /n n ^ V 

\, tim ^\--^.«.[- lf -)) + ^>,s,n^— Jj (io> 

m = 0,1,2, ... ,— 

[0031] Note that the pitch of the teeth around the encoder tone wheel 

circumference is not constant because of variations in manufacturing tolerances, 
and therefore unbalanced variations for every revolution of tire occurs. Frequencies 
associated with these variations are referred to as pole pitch frequencies, and a pole 
pitch error is defined as the maximum tolerance on the teeth of the encoder tone 
wheel. 

[0032] As the tire rotates, the pitch error will periodically create additional 

vibration in the frequency spectrum. These additional frequencies make it more 
difficult to find a peak torsional vibration in the angular domain. Therefore, in 
accordance with the invention, the pole pitch frequencies of the single spectrum are 
eliminated before proceeding with the detection algorithm or process. 
[0033] Specifically, in the angular domain, pole pitch vibration occurs where 

the angular frequency n is equal to a multiple of cycles I 2k radian . Although there is 
no relation between pole pitch vibration and the vehicle speed in the angular 
domain, it is easier to eliminate pole pitch vibration in the angular domain than in the 
time domain because the vibration occurs at fixed angular frequencies in the angular 
domain. 

[0034] FIG. 3 shows an angular frequency spectrum with and without pole 

pitch frequencies as a vehicle is driven on a regular road at the speed of about 80.0 
kph, illustrating that pole pitch error distorts the frequency spectrum. 
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[0035] Although the spectrum in FIG. 3 is shown in the angular frequency 

domain, for vehicle vibration analysis, it is more convenient to express the amplitude 
of surface profiles in terms of the time frequency in Hz rather than in terms of the 
angular frequency in cycles I 2 n radian since the vehicle vibration is a function of 
time. Hence, the angular frequency spectrum n in cycles/ ix is transferred to the 
time frequency spectrum in Hz with the following expression: 

f Hz - CI {cycles I2tt radian )xa> (2n radian I sec) (11) 

where 

n - Angular frequency (cycles 1 2n radian ) 

co - Angular velocity of a wheel {2n radian I sec) . 
And the amplitude of the angular profile in terms of the angular frequency a 6 (Q) is 
transferred to the temporal frequency^ ,(/) , according to the expression: 

= ^ (12) 

CO 

where 

A t = amplitude of the temporal frequency spectrum, and 

A Q = amplitude of the angular frequency spectrum. 
[0036] FIG. 4 shows in the time frequency domain that after the pole pitch 

errors are eliminated, the peak for-aft torsional vibration shifts from 47.0 Hz to 35.0 
Hz for a tire pressure loss of 1.2 bars for a vehicle driven over a regular road along a 
substantially straight line at a speed of about 80.0 kph. FIG. 5 shows the peak for- 
aft torsional vibration shifts from 42.0 Hz to 35.0 Hz for a tire pressure loss of 1.2 bar 
for a vehicle driven over a rough road at a speed of about 40 kph. 

10 



Attorney Docket No. 10543-074 



[0037] As can be seen FIGs. 4 and 5, even after eliminating pole pitch error, a 

single frequency spectrum still includes noise or deviation. Therefore, to eliminate 
deviation of the single spectrum, an averaging spectrum s p (mACl) is introduced by 

combining the continuous K single spectra according to the expression: 



s p (mAQ.) = ]T s k (mAQ ) m = 1,2,3, M - 1 



(13) 



*=0 



[0038] Curve-fitting is then applied to the series of averaged frequency 

spectra. Specifically, a peak torsional vibration frequency F 0 (cycleslln) in the 
angular domain is curve fitted according to the expression: 

f n = — — (14) 



where, 



2 x A, 



AY, -I 


2>„ 3 , 




m —M | 


m=A/, 


m=AY , 


M 2 -\ 


M 2 -\ 


AY 2 -1 








m=M i 


£ y m 

m—M | 


A/ 2 -l. 




m =A/ 1 


M 2 -M 



(15a) 



m s m 



M 2 -\ 
m=A/, 

_ A* 2 -1 

m-M K m=M t m=M , 

AY, -I A/, -I 

m-M | w=AY, 



ivi — i yvv -> — 

m=M | m=M t 
AY,-1 M 2 -\ 



M 2 - A/, 



(15b) 



and 

jc„, = mAQ 

(16) 

J/, =S„(mAQ) m = M,,M,+l, ...,M 2 -1 
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[0039] The range of the curve-fitting is determined by the values of M y and 

M 2 , which are functions of the angular velocity co . 

[0040] A 'BAD' spectrum is rejected if the fitted curve regresses to a straight 

line (i.e., if A a =0 ), if the peak position is out of the pre-defined range, or if the 
calculation overflows. 

[0041] As mentioned above, it is more convenient to characterize vehicle 

vibrations in the time frequency domain rather than in the angular frequency domain. 
Therefore, after calculating the peak torsional vibration from the frequency spectrum, 
the process maps the peak torsional vibration from angular domain to the time 
domain. That is, the process continues to calculate a series of frequency 
observations and determines an initial torsional frequency (f 0 ) in the time frequency 

domain based on the observations. Specifically, the peak torsional vibration is 
mapped to the frequency domain (Hz) according to: 

fo=F 0 xa> 0 (17a) 

For example, assuming & = in x / and the tic time is ju = 4 x 10 ~ 6 sec , 
then the angular velocity is 

(17b) 



(Z A/„)//*4xl0 



(2/r / sec) 



where / is a predefined number of sampling cycles. 

[0042] FIG. 6 illustrates the results of curve-fitting applied to the frequency 

spectrum generated for a vehicle driven over a rough road at the speed of about 
60.0 kph after mapping to the time frequency domain. As shown, the peak 
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frequency shifts from 40.0 Hz to 35.0 Hz during a tire pressure loss from 2.7 bars to 
1.5 bars. 

[0043] After calculating the initial frequency, the process makes adjustments 

to the frequency over time. Even if the vehicle changes speed, the process 
continues to calculate the single spectrum, the averaging spectrum, and the peak of 
angular frequency in the angular domain and maps the peak F j in the angular 
domain to the time domain /, ( / = 1,2,3, ... ) by the aforementioned process, that is, 

f k =F, *G) n i = l,2,3,...,/,... (18) 
[0044] A filtering calculation is also applied to the series of observed values. 
Moreover, when / reaches the predefined value / , the process outputs and 
continues to update the peak frequency of the torsional vibration in real time. 
[0045] The aforementioned process is able to detect pressure losses in more 

than one tire, since, for example, tire pressure loss in the front left tire causes 
vibration shifting in that tire but does not impact the performance of the front right tire 
and vice versa. Similarly, vibration of the rear left tire is isolated from the vibration of 
the rear right tire and vice versa. 

[0046] As an example, the process was implemented in a front-wheel drive 

vehicle driven on a rough road at about 60 kph, in which the tire pressure of the front 
right tire was 2.3 bars and that of the rear right tire was 2.1 bars. The tire pressure 
of the front left and rear left tires were changed to 1.5 bars, 1.9 bars, and 2.5 bars. 
As shown In FIG. 7, pressure loss in the front left tire caused the peak for-aft 
vibration to shift from 47.3 Hz to 40.9 Hz. As shown in FIG. 8, pressure loss in the 
rear left tire shifted the peak for-aft vibration from 47.0 Hz to 43.8 Hz. 
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[0047] Note also in FIGs. 7 and 8 that the amplitude of the peak vibration 

associated with the rear left tire is less than that of the front left tire. This is 
expected, since in a front wheel drive vehicle torsional compliance is provided by the 
front wheel axles and other powertrain components. The non-driven rear wheels, on 
the other hand, have very little connected components and therefore there is no 
significant coupled torsional compliance associated with the rear wheels. 
[0048] FIG. 9 illustrates that there is no significant shift in the peak for-aft 

vibration associated with the front right tire. That is, the front left axle vibration does 
not impact the performance of the front right axle. Similarly, FIG. 10 shows that 
there is no significant shift for the rear right tire. 

[0049] The above described process is illustrated in FIG. 11 as a state 

diagram 100 and in FIG. 12 as a flow diagram of a process of steps 200. Referring 
in particular to FIG. 12, after initiating in step 202, the process 200 determines the 
discrete (Dt) contribution in step 204 from which the real and imaginary components 
are obtained to calculate the spectrum in step 206. 

[0050] In step 208, after creating the single spectrum, the process 200 

eliminates the pole pitch error and in step 210 combines the continuous single 
spectra to an averaging spectrum in the angular domain. In step 212, the process 
200 curve fits the averaged spectrum and maps the peak frequency of the curve- 
fitted spectrum from the angular domain to the time domain in step 214. 
Subsequently, the process 200 filters the peak frequency and makes long-term 
adjustments to the peak frequency over time in step 216. As mentioned above, the 
process 200 continues to calculate the torsional vibration frequency, even when the 
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vehicle speed changes, and applies statistical analysis to the series of observed 
data to extract the best value of the estimated frequency. This update is used to 
detect the shift in the torsional vibration frequency, which is related to the tire 
pressure loss of one or more tires in step 218. The results are displayed in step 
220, and the analysis subsequently ends. 

[0051] Other embodiments are within the scope of the following claims. 
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